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ABSTRACT 


The  effects  of  a  sloping  bottom  on  acoustic  transmissions,  between  a  source 
and  receiver  at  arbitrary  but  fixed  locations,  are  investigated  using  ray 
theory.  An  isospeed  channel  is  assumed,  and  bottom  angles  up  to  about  3°  are 
considered-  Sloping  bottom  influence  on  per-ray  quantities,  including  travel 
time  and  transmission  loss,  are  examined  for  cw  transmissions .  Significant 
variations  are  shown  to  occur,  such  as  travel  time  changes  of  more  than  200  ms 
over  ranges  of  about  six  km.  Per-ray  transmission  loss  is  found  to  be 
influenced  strongly  by  bottom  slope,  the  amount  of  influence  depending  upon 
source-receiver  bearing  and  the  bottom  loss  model  used.  Variations  of  more 
than  20  dB  are  demonstrated.  Effects  of  a  sloping  bottom  on  the  total 
acoustic  field  are  examined  also,  and  the  results  compared  with  those  for  a 
horizontal  bottom.  Finalxy,  a  simple  model  of  a  shallow  water  front  is 
superpcsed  over  the  sloping  bottom,  and  travel  time  is  investigated.  The 
sloping  bottom  effect  can  induce  travel  time  changes  more  than  300%  larger 
than  the  frontal  effect  for  different  source-receiver  geometries  and  be atom 


inclinations • 
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INTRODUCTION 

Shallow  water  in  the  world's  oceans  occurs  along  continental  margins  and 
in  many  of  the  seas,1  Continental  margins  are  subdivided  into  the  continental 
shelf,  a  region  of  gently  sloping  bottoms  with  slopes  of  0°  to  3°,  and  the 
continental  slope  or  shelf  break  region,  where  water  depth  increases  rapidly 
and  bottom  inclination  can  be  as  much  as  Acoustic  transmission  in 

shallow  water  is  signif icantly  different  from  deep  water  transmission  with  the 
former  being  characterized  by  stronger  bottom  interaction  and  somewhat  shorter 
ranges  (see,  for  example,  Ref.  3).  Acoustic  models  of  the  shallow  ocean  vary 
in  complexity  from  a  simple  constant  depth,  isospeed  medium  to  one  in  which 
depth,  sound  speed,  and  bottom  properties  vary  spatially. ^ 

A  variety  of  papers  have  shown  certain  consequences  of  shallow  ocean 
environmental  features  on  acoustic  transmissions.  For  example,  both  ray  and 
mode  theories  were  used^'8  to  indicate  influences  of  sloping  bottoms  on 
signal  intensity  and  on  horizontal  trajectories  of  propagating  modes.  A 
method7  for  studying  the  horizontal  propagation  of  local  normal  modes  has  been 
applied  to  a  wedge  shaped  duct,  which  affects  individual  mode  trajectories 
and  eigenvalues  differently.  Accurate  analytic  approximations  from  ray 
theory,  for  ray  geometry,  tra^ai  time,  and  spreading  less,  have  been 
obtained,8' 9  but  only  for  the  strictly  two-dimensional  propagation  directly 
up-or  downslope.  Transmission  losses  and  shadow  zones  for  an  isospeed  wedge 
shaped  duct  have  been  described  for  medium  range  propagation  using  geometric 
acoustics  and  normal  modes.10  The  three-dimensional  nature  of  acoustic 
propagation  induced  by  variable  bottom  tocography  has  been  shown  to  strongly 
influence  horizontal  ray  trajectories  and  shadow  zone  locations.11'1-  Recent 
studies  have  concentrated  on  application  of  the  parabolic  approximation  to 
shallow  water  sloping  bottom  oceans,1^1  4  to  show  modal  cutoff  for  ups lope 
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propagation . 1  5  Also,  modal  coupling  has  been  considered , 1  6  an<3  different  methods 
for  wave  propagation  over  an  inclined  bottom  have  been  developed1'  and 
adapted . 1 ® ' 1 ^  Relationships  among  propagation  theories  have  been  clarified  in 
detail  (as  In 'Ref.  20).  In  addition,  some  deep  water  studies  have  used  ray  theory 
to  demonstrate  inclined  bottom  effects  between  fixed  sources  and  receivers . r 22 

The  purpose  of  this  paper  is  to  study  sound  transmission  in  shallow  water 
over  a  sloping  bottom  when  source  and  receiver  are  located  at  arbitrary 
positions  in  the  sound  channel,  but  at  relatively  short  range.  We  consider  cw 
signals  of  sufficiently  high  frequency  to  permit  the  use  of  ray  theory.  The 
influence  of  the  sloping  bottom  on  acoustic  quantities  is  developed  in  a 
constant  sound  speed  medium  with  a  Mackenzie  bottom  model. ^3  other  bottom 
models  may  be  used  in. shallow  water  (see,  for  instance,  Ref.  24)  with  similar 
results.  A  sound  channel  with  horizontal  surface  and  bottom  dictates 
two-dimensional  ray  paths  in  the  vertical  plane  containing  the  source  and 
receiver.  Imposition  cf  a  plane,  sloping  bottom  necessitates  consideration  of 
three-dimensional  ray  paths  for  other  than  directly  up-  or  downslope,  as  is 
well  known. 

The  principal  novelty  of  this  paper  is  the  development  of  analytic  approx¬ 
imations  for  three-dimensional  ray  geometry,  travel  time,  and  per-ray  and 
total-field  loss  between  a  fixed  source  and  receiver.  Our  study  of  ray  geometric 
properties  and  of  certain  acoustical  effects  of  three-dimensional  ray  paths  pro-= 
duces  new  results.  For  instance,  our  formulas  facilitate  interpretation  of 
acoustic  variations  which  result  from  changes  in  source-receiver  location  and 
bottom  inclination.  Another  feature  of  this  paper  is  an  appraisal  of  the  rela¬ 
tive  acoustic  significance  of  a  sloping  bottom  and  a  shallow  ocean  front,  using 
results  of  another  study. In  numerical  examples  we  use  frontal  parameters 
which  are  similar  to  those  of  the  Slope  Front  in  the  North  Atlantic  Ocean.  Enough 
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parameter  variations  are  considered  so  that  our  results  are  applicable  to  a 
variety  of  ocean  environments.  Finally,  we  remark  that  detailed  comparisons 
of  our  results  with  those  in  the  previous  studies  cited  above  are  not  attempted 
in  this  paper*  This  is  because  different  acoustical  quantities  are  typically 
emphasized,  and  because  our  results  are  applicable  to  shorter  range  trans¬ 
missions,  while  previous  studies  are  primarily  applicable  to  longer  ranges. 
However,  it  should  be  noted  that  there  is  strong  qualitative  agreement  between 
level  curves  of  transmission  loss  in  this  paper  and  those  calculated  by 
Bradley1^,  which  were  obtained  using  normal  mode  theory  and  which  apply  at 
longer  ranges. 

In  Sec.  I  we  describe  ray  geometry  determined  by  regular  perturbations 
from  an  isospeed  sound  channel  with  horizontal  surface  and  bottom.  It  is  shown 
how  ray  geometry  depends  upon  bottom  inclination,  source-receiver  bearing 
relative  to  the  maximum  gradient  of  bottom  slope  and  range.  Results  are  used 
in  Sec.  II  to  express  travel  time  and  spreading  loss  in  terms  of  bottom  angle 
and  bearings.  Then,  a  Mackenzie  bottom  model  is  used  in  an  examination  of 
per-ray  transmission  loss.  Incoherent  total-field  transmission  loss  is 
considered  also,  and  the  sloping  bottom  influence  is  discussed.  Section  III 
studies  the  effect  of  a  sloping  bottom  superimposed  on  propagation  through  a 
simple  frontal  model.  Ray  geometry  and  travel  time  expressions  are  derived, 
and  the  relative  significance  of  the  sloping  bottom  and  front  are  illustrated. 
Finally,  we  summarize  our  results  in  Sec.  IV. 
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I.  MODEL  DESCRIPTION  AND  RAY  GEOMETRY 

The  bottoms  on  the  continental  shelf  and  in  shelf  break  regions  are 
commonly  and  most  easily  modelled  as  planes  with  angles  of  inclination  of  up 
to  about  5Q.4.,8,9  In  this  paper,  we  model  the  ocean  above  the  continental 
shelf  as  a  channel  bounded  above  by  a  horizontal  surface  and  below  by  a  plane 
having  an  angle  of  inclination  a.  Subsequently,  we  will  use  radian  measure 
for  this  and  other  angular  quantities.  In  order  to  concentrate  on  the 
sloping  bottom  effect,  we  initially  suppress  oceanic  variations,  taking  the 
water  to  be  stationary  and  to  have  constant  sound  speed  c.  The  influence  of 
the  sloping  bottom  will  be  investigated  using  ray  theory  and  a  Mackenzie 
bottom  interaction  model. 23 

Suppose  that  an  omnidirectional  sound  source  S  and  a  point  receiver  R  are 
located  at  arbitrary  depths  and  hR  in  an  ocean  channel  of  depth  d^  at  R. 

We  establish  a  right-handed  Cartesian  coordinate  system  with  origin  on  the 
surface  over  R,  with  depth  z  increasing  downward,  and  with  x  increasing 
directly  upslope  (shoreward),  as  shown  in  Fig.  1(a).  The  bottom  angle  a  is 
nonnegative  and  measured  from  the  horizontal,  and  the  relative  bearing  0  of  5 
from  R  is  measured  positively  clockwise  from  the  y-axis'bas  shown  in  Fig.  1{b). 
It  is  the  angle  0(-tt  <  0  <_  tt  ) ,  then,  which  specifies  propagation  in  a  conver¬ 
gent  or  divergent  channel.  The  range  R  is  the  horizontal  distance  separating 
S  and  R.  We  use  N  to  designate  the  total  number  of  bottom  reflections  between 
S  and  R  for  a  given  ray,  and  we  number  bottom  reflections  sequentially  from  5 
to  R  by  1,2,...,N.  A  ray  lobe  is  that  portion  of  a  ray  path  between 
successive  bottom  reflections,  and  the  lobe  is  numbered  with  the  smaller  of 
its  bottom  reflection  numbers  (i.e.,  lobe  k  is  than  portion  of  the  ray  path 
between  bottom  reflections  k  and  k  +  1).  Lobe  zero  is  the  (typically  partial) 
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lobe  between  S  and  the  first  bottom  reflection,  while  lobe  N  is  the  (typically 
partial)  lobe  between  the  last  (Nth)  bottom  reflection  and  R.  We  designate 
the  water  depth  at  S  by  dg  and  the  depth  at  bottom  reflection  k(1  <  k  <  N)  bv 
hkN  • 

The  altitudinal  angles  of  a  ray  at  S  and  R  are  0SN  and  0RN,  respectively, 
measured  positively  downward  from  the  horizontal.  In  lobe  k  (1  <  k  <  N-l  )  , 

©kN  is  the  positive  altitudinal  angle  of  the  ray  after  its  surface  reflection; 
©oN  is  the  altitudinal  ray  angle  in  lobe  zero  as  the  ray  path  approaches  the 
first  bottom  reflection;  0NN  is  that  angle  in  lobe  N  after  surface  reflection, 
if  such  a  reflection  occurs.  For  arbitrary  source  and  receiver  depths  and 
fixed  N,  there  will,  in  general,  be  four  rays  connecting  S  to  R.  We  introduce 

the  parameter  as(aR)  which  is  +1  if  0SN  (0RN)  >  0  and  is  -1  if  0SN  (0^)  <  0. 

We  relate  altitudinal  angles  by 

0SN  =  aS0oN  (la) 

and 

0RN  =  aR0NN  (1b) 

In  Fig.  1,  we  show  schematically  a  ray  with  N  =  3,  for  which  as  =  -1  and  oR  = 

+  1  .  Each  ray  lobe  lies  in  a  single  vertical  plans  and  specular  reflection 
insures  that,  in  general,  adjacent  lobes  lie  in  different  vertical  planes, 
giving  rise  to  three-dimensional  ray  paths.  The  azimuthal  angle  in  lobe  k(0  < 

k  _<  N)  of  a  ray  is  measured  in  the  same  way  as  0.  The  projection  ci  tie 

ray  path  in  lobe  k  onto  a  horizontal  plane  is  called  LkN.  In  general, 
N-dependent  quantities,  such  as  8SN  and  hkN,  are  functions  of  os  and  aR.  Fcr 
simplicity,  however,  we  will  suppress  this  dependence  in  subsequent  notation. 
Further,  nondlmensional  quantities  will  be  designated  by  a  caret  ( - )  ,  while 
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quantities  associated  with  a  sound  channel  having  a  horizontal  bottom  of 

depth  dR  will  be  indicated  by  an  overbar. 

For  fixed  source  and  receiver  locations,  fixed  N,  and  specified  as  and 

oR,  we  facilitate  the  determination  of  the  geometry  of  a  ray  by  specifying  the 

unit  vector  tangent  to  the  ray  at  5,  Ag^,  and  at  R,  ARN.  In  lobe  k(  0  <_  k  <_  N)  , 

the  unit  vector  tangent  to  the  ray  approaching  the  bottom  reflection  is  A*  . 

“kN 

Although  A*  is  not  actually  present  on  the  physical  ray  when  aR  =  -1  ,  it  is  a 
computational  convenience;  a  similar  statement  applies  to  AoN  when  as  =  +1  . 
Unit  ray  tangents  may  be  written  as 


and 


cos  0SN  sin  4»oN' 
-cos  0SN  cos  <j)oN 
sin  0SN 


(2a) 


4rn 


cos 


yRN 


sin  (j) 


-cos  0R^ 
sin  0RN 


NN 

^NN 


(2b) 


(^kN 


cos  0kN  sin. 

-cos  0kN  cos  <}>kN 
-sin  0kN 


0  <  k  <  N 


(2c) 


The  effect  of  a  boundary  reflection  on  the  unit  tangent  vector  to  a  ray 
is  normally  written  as  a  vector  equation,  with  the  reflected  tangent  vector 
expressed  as  a  linear  combination  of  the  incident  unit  tangent  vector  and  a 
vector  normal  to  the  boundary. This  linear  combination  of  vectors  is  a 
linear  transformation  of  the  incident  ray  tangent,  and  can  be  represented  by 
a  matrix.27  The  matrix  of  transformation  T.j  at  a  surface  reflection,  and 
that  at  a  bottom  reflection,  Tc,  can  be  written  as 


i '  ° 

It  -  °  ' 

\  0  o 


(3a) 


and 
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(cos  2a  0  -sin  2a\ 

0  1  °  J  .  (3b) 

-sin  2a  0  -cos  2a/ 


Ray  tangents  can  now  be  related  by 


A'=tA  ,  0  <  k  <  N  ,  (3c) 

~kN  ~T  ~kN  —  ~ 


and 


A'  =  T  T  A'  ,  0  <  k  <  N-1  .  (3d) 

(k+1 )N  T  3  kN 


Thus,  we  can  relate  A'  to  a  ray  tangent  near  R  by 

~kN 


A'  =  (t  T  )N-k  A'  ,  0  <  k  <  N-1  .  (4) 

~kN  ~B~T  ~NN  —  — 


The  ray  tangent  at  S  (or  R) ,  ASN  (or  ARN) ,  is  AqN  (or  A^)  if  as  (or  crR)  =  -1  . 


It  is  necessary  to  determine  ray  geometry  ir  order  to  calculate  travel 


time,  spreading  loss,  and  bottom  loss.  We  first  seek  to  determine  the 


altitudinal  and  azimuthal  angles  of  a  ray  tangent  at  S  and  R,and  these 


quantities  are  than  specified  at  any  point  on  the  ray  path  through  Eq„  (4). 


We  write  the  x  and  y  components  of  range  in  the  form 

N 

R  sin  0  =  l  sin  ^  (5a) 

0 

ana 

N 

R  cos  8  =  1  L;<n  cos  <t>kN  .  ,  (5b) 

0 

Now,  the  are  related  to  water  depth,  ray  angles,  and  bottom  inclination 


by  the  equations 
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LoN  =  (ds-ashs)  ( tan  0oN+tan  a  sin  <{>oN)  1  ,  (6a) 

LkN  =  2hkN(tan  0kN+tan  a  sin  i|>kN)-1  ,  1  _<  k  _<  N-1  ,  (6b) 


and 


lNN  =  (dR+aRhR)  (tan0NN-tan  a  sin  <j>NN)  1  .  (6c) 

Further,  water  depths  at  successive  bottom  reflections  can  be  determined 
recursively  from 

h(k+1)N  32  hkN(  tanGkN-tan  a  sin  <J>kN) 

x  (tan  0kN-i-tan  a  sin  <j>kN)-l  ,  1  _<  k  <_  N-1  ,  (6a) 

with  h1N  being  given  by 

h1N  =  (ds  tan  0qN  +  ashs  tan  a  sin  <j>oN) 

x  (tan  0oN+tana  sin  •  (6e) 

When  Eos.  (6)  are  substituted  into  Eqs .  (5)  and  (4).  we  have  a  system  of  non¬ 
linear  coupled  equations  in  the  four  angles  0oN,  ^oN'  4>nn*  Equations 

(1)  then  give  0SN  and  0rn»  Since  the  system  cannot  be  solved  analytically,  we 
proceed  to  perturb  off  known  results^  for  an  isospeed  channel  with  constant 
depth  dR  (i*e.,  a=Q)  .  For  a  <<  1 ,  we  write 

®kN  =  "**  ^0kN  /  (7a) 


and 


<t>kN  =  S  +  A<frkN  ,  0  £  k  <_  N  ,  (7b) 

where  0^  is  known  from  Ref.  3  to  be 


9 


tan  0^  =  [  2N  -  ashs  •+  aR^R  +  2H(hg-1)N  tan  ajH-"1  .  (8) 

In  Eq.  (8),  depth  and  range  parameters  with  carets  are  nondimensionalized  with 

respect  to  d  ,  h  =  h  d”1  ,  h  =  h  d-1  ,  and  R  =  Rd_1  ,  and  H  is  the  Heaviside 

RSSRRRR  R 

function.  The  Heaviside  function  permits  inclusion  of  the  last  term  in  Eq. 

(8)  whenever  source  depth  exceeds  dR. 

When  we  substitute  Eqs .  (7)  into  Eqs .  (2)  and  (4)  and  linearize  in  A0kN, 
A<t>kN  and  a'  we  obtain  approximations  for  A0kN  and  A(J>kN  in  terms  of  A0oN  and 
A<j>oN.  The  approximations  for  A0kN  and  A<j>kN  are  substituted  first  in  Eqs.  (6), 
and  then  in  Eqs.  (5),  with  linearization  in  small  quantities  at  each  step. 
Eventually  we  are  able  to  approximate  angular  changes  in  the  four  ray  angles 
at  S  and  R,  in  the  form 

A®oN  “  -2Na(N+aRhR)R-1  sin  3  cos  0N  sin  0N  ,  (9a) 

AOjnj  =  -2Na[ (N+aRhR)R-1  sin  3  cos  0N  sin  0N  -  2N  sin  3j  ,  (9b) 

A<{>oN  *  2Na(N+aRhR)R_1  cos  3  ,  (9c) 

and 

A^njj  =  2Na[  (N-!-cRhR)R-1-tan  0N]cos..  3  .  (9d) 

Equations  (9)  are  accurate  provided  0N  is  not  too  close  to  0  or  tt/2  rad, 
so  that  unperturbed  rays  interacting  with  the  bottom  must  not  be  too  shallow 
or  too  steep.  In  order  to  assess  the  accuracy ,  approximations  from  Eqs.  (9) 
were  compared  with  numerical  solution  of  the  system  of  equations.28  For 
example ,  with  surfaced  S  and  Rf  R  =  5  km,  a  =  0.5°,  8  =  -45°,  dR  =  300  m,  and 
rays  with  N  =  1,2,... ,5,  approximations  to  altitudinal  angles  showed  relative 
errors  of  less  than  4%,  while  bearing  approximations  were  accurate  to  within  6%. 
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The  altitudinal  angles  of  perturbed  rays  transmitted  to  a  receiver  ups  lope 
( 0  <  0  <  rr )  are  shallower  than  their  unperturbed  ray  counterparts,  i.e.  A0qN  < 
0;  yet  the  altitudinal  arrival  angles  at  R  are  steeper  because  of  bottom 
interaction.  Similarly,  rays  launched  downslope  (-7T  <  0  <  0)  start  more 
steeply  than  unperturbed  rays,  A0qN  >  0,  and  arrive  at  shallower  angles.  When 
S  is  directly  across-slope  from  R  (0  ~  0,70,  altitudinal  angles  are 
insignificantly  affected  by  bottom  inclination  as  shown  in  Eqs .  (9a)  and  (9b). 
The  azimuthal  changes  indicated  in  Eqs.  (9c)  and  (9d)  show  that  the  tendency  of 
the  sloping  bottom  to  reflect  a  ray  slightly  downslope  with  each  reflection  is 
countered  with  a  launch  azimuth  more  upslope  than  the  source  bearing  3.  Thus, 
we  see  that  A4>oN  is  positive  for  S  in  the  first  and  second  quadrants  (y  >  0)  , 
and  negative  for  5  in  the  third  and  fourth  quadrants,  of  the  xy-plane  shown  in 
Fig.  1.  Of  course,  there  is  no  azimuthal  variation  when  the  source  is  directly 
upslope  (0  =  tt/2)  or  directly  downslope  (0  =  tt/2)  . 

II.  TRAVEL  TIME  AND  TRANSMISSION  LOSS 

First,  we  develop  and  examine  an  expression  for  the  travel  time  for  a 
ray  with  N  bottom  reflections  (and  with  a  specification  of  as  and  aR)  .  If  SN 
is  the  arc  length  of  the  ray,  then 

N 

TN  =  SN/c  =  c-1  l  LkN  sec  Qj^  .  (10a) 

0 

If  we  define  Atn  to  be  the  change  in  travel  time  from  that  in  an  ogean  with 
horizontal  bottom  of  depth  dR,  T^,  then 

ATn  =  -  T^j  ,  (10b) 

where 

TN  =  (R  sec  0^)c-1  . 


(10c) 
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Using  the  angle  changes  of  Eqs .  (9)  to  approximate  the  L^N  in  Eqs .  (5),  and 

substituting  into  Eqs.  (10a)  and  (10b),  it  can  be  shown  that,  to  first  degree 
terms  in  a, 

A-Tn  =  2NaTN{  sin  0 [ N+aRhR] r[ R2+( 2N  +  ashs  +  Cr^r)2]”1}  •  (  1  0d ) 

A  suggestion  of  the  accuracy  of  our  approximation  can  be  given  for  the  para¬ 
meter  values  used  in  Sec.  I.  In  this  case,  the  approximate  travel  time 
calculated  from  Eqs.  (10b)  and  ( 1 Od )  has  been  shown  to  vary  with  N  from  3.1% 
to  11.1%  of  that  predicted  by  Eq.  (10a),  in  which  results  of  the  aforementioned 
numerical  solution  for  ray  geometry  are  used. 

To  illustrate  travel  time  variations  here,  and  variations  in  other 
acoustical  quantities  subsequently,  we  shall  use  the  parameter  values  dR  =  300 
m,  c  =  1500  ms-1,  and  a  =  2°  .  For  convenience,  we  shall  take  S  and  R  in  our 
numerical  examples  to  be  located  on  the  surface.  Other  source  and  receiver 
depths  have  been  shown  to  give  rise  to  similar  behavior  in  acoustical 
quantities.  We  note  that  S  and  R  on  the  surface  correspond  to  the  values  c3  = 

+  1  and  oR  =  —  1  • 

Figure  2  displays  level  curves  of  T5  from  Eqs.  (10b)-(1Cd)  for  R  fixed  at 
the  origin  and  for  5  at  (x,y).  The  level  curves  are  clearly  skewed  in  the 
ups  lope  direction.  The  observed  x-axic  symmetry  can  be  seen  to  come  from  the 
sin  0  term  in  Eq.  (lOd).  Also,  the  near-circularity  of  the  level  curves  is 
explained  by  considering  T5  from  Eq0  (10b)  as  the  polar  equation  of  a  Limaccr 
of  Pascal,  r  =  b  +  a  sin  0,  with  a  r<  b  after  substituting  Eqs.  (10c)  and 
(lOd).  Travel  times  for  sources  equidistant  from  R  are  generally  longer  when 
propagating  upslope  (0  >  0)  than  when  propagating  downslope,  with  longest 
travel  time  corresponding  to  directly  upslope  transmission  (0  -  tt/2)  .  When  S 
is  located  in- water  of  the  same  depth  as  that  at  R  (0  =  0  or  it),  travel  times 
from  sources  equidistant  from  R  are,  of  course,  equal. 


In  Fig.  3,  we  allow  5  to  assume  successive  positions  along  each  of  four 
paths ,  P1,.../P4,  which  are  shown  in  Fig.  2  and  which  have  the  same  closest 
point  of  approach  (CPA)  to  the  source,  5  km.  Path  Pj  represents  a  track 
parallel  to  the  shore  upslope  from  R  (sound  propagates  downs lope ) .  Path  ?2 
is  similar,  except  that  the  track  is  downslope  from  R  and  the  signal  propa¬ 
gates  upslope.  Path  P3  represents  a  track  directly  upslope.  A  track  at  45° 
across  the  slope  is  represented  by  path  P4.  The  symmetry  of  travel  time  T5 
about  the  x-axis  can  be  seen  on  paths  P-j  and  P2  in  Fig.  3.  In  contrast,  the 
asymmetry  of  T5  about  the  y-axis  on  path  P3  is  seen  from  the  shift  of  the 
minimum  away  from  CPA.  This  is  caused  by  the  fact  that  arc  length  is  a 
minimum  for  a  source  location  slightly  upslope  from  3=0.  When  compared  to  a 
constant  depth  ocean,  a  =  0,  it  can  be  seen  that  a  source  on  a  track  upslope 
of  R ,  such  as  P-j  ,  generates  shorter  travel  times  than  one  on  a  path  downslope 
from  R,  such  as  P 2,  by  hundreds  of  ms.  On  path  P3,  T5  is  larger  than  T5  (a=Q) 
initially,  but  intersects  T5  at  CPA  and  becomes  less  than  the  unperturbed 
travel  time  as  S  moves  upslope  from  R.  On  the  diagonal  path,  P4,  travel  times 
are  greater  than  all  other  paths  initially,  become  less  than  those  on  path  P3 
after  intersection  with  that  path,  and  then  approach  the- a  =  0  result  as  S 
approaches  the  y-axis.  We  note  that  on  the  y-axis  both  T5  and  T5  should  be 
equal. 

The  variation  of  travel  time  from  the  horizontal  bottom  state  is  more 
easily  seen  in  Fig.  4  where  we  present  level  curves  of  the  change  in  travel 
time,  AT5,  given  by  Eq.  (IQd).  We  observe  that  AT5  >  0,  or  T5  >  T5 ,  when  S  is 
downslope  from  R  (upslope  propagation),  as  expected.  Besides  the  casual 
observation  that  | AT5  j  increases  with  range,  we  see  clearly  from  Fig.  4  that 
AT5  changes  most  rapidly  when  progressing  directly  upslope  along  the  x-axis. 
Near  R,  AT5  changes  as  much  as  50  ms  per  km  when  propagating  upslope.  Thus, 
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bottom  inclination  strongly  influences  travel  time,  the  magnitude  of  this 
influence  being  determined  by  the  relative  bearing  of  5  from  R. 

Geometric  spreading  loss  for  a  ray  with  N  bottom  reflections  and  for 

ag  and  aR  specified,  can  be  derived  from  the  equation^ 


Sln  9RN  (  3x  3y  _  8x  3y  j 
cos  ®SN  9®SN  3<^SN  ^SN  30SN 


(11a) 


where  the  partial  derivatives  are  evaluated  at  x  =  y  =  0,  and  where 


and 


N 

x  =  R  sin  0  -  l  LkN  sin  <j>kN 
0 

N 

y  =  R  cos  0  -  l  LkN  cos  <j>kj}  . 
0 


( 11b) 


(11c) 


Performing  the  operations  indicated  in  Eq.  (11a)  on  Eqs .  (11b)  and  (11c),  and 
making  substitutions  from  Eqs.  (1)-(5)  and  (9),  we  find  that,  to  lowest  order 
terms  in  a,  spreading  loss  is  approximately  the  square  of  arc  length: 


L  =  s2  =  [cT  ]2  .  (lid) 

N  N  N 

This  equation  shows  that  spreading  loss  varies  as  the  square  of  travel  time. 
Therefore,  at  fixed  R,  bottom  slope  causes  LN  to  be  larger  for  signals 
transmitted  ups lope  than  for  signals  propagated  downs lope,  as  was  seen  for 
travel  time.  Although  both  bottom  inclination  and  bearing  of  S  from  R  affect 
Ln,  the  variations  from  the  horizontal  bottom  result  are  less  than  2  dB  for 
typical  bottom  parameters. 

In  contrast  to  spreading  loss,  bottom  loss  is  known  to  be  particularly 
significant  in  shallow  water  propagation.24  Thus,  bottom  slope  can  be 
expected  to  have  an  impact  on  transmission  loss,  -20  log^Q  AN,  where 

N 

a  =  n  a  l  -1/2 

N  !  kN  N 


(12) 
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is  amplitude  and  8^  is  the  reflection  coefficient  at  bottom  bounce  k  of  a  ray 
with  N  bottom  reflections  (and  <Jg,  aR  specified)  .  We  calculate  using  a 

Mackenzie  bottom  model  and  our  previously  developed  geometric  approximations 
for  0^  and  •  For  numerical  examples,  the  parameters  used  are  typical  of  a 
sand  bottom,  so  that  we  take  bo ttom-to-water  sound  speed  ratio  to  be  1.1653, 
bottom- to-wa ter  density  ratio  to  be  1.9522,  and  an  attenuation  parameter  of 
0.01 58. 23  Other  bottom  parameter  values  have  been  shown  to  generate  similar 
acoustical  results. 

Level  curves  of  per-ray  transmission  loss  appear  in  Fig.  5,  again  for 
surfaced  source  and  receiver  and  N  =  5.  We  note  that  the  loss  drops  as  S-R 
range  increases  to  about  6  km,  and  then  increases  with  R,  although  spreading 
loss  increases  monotonicaily  with  R.  This  transmission  loss  behavior  is 
attributed  to  the  bottom  loss  model,  in  which  the  steep  angles  of  incidence  at 
relatively  short  ranges  cause  much  greater  bottom  loss  than  that  due  to 
spreading.  As  range  increases,  however,  the  angle  of  incidence  of  a  ray 
decreases,  causing  bottom  loss  to  decrease  while  spreading  loss  increases. 

The  transmission  loss  ultimately  increases  with  range,  when  spreading  loss 
becomes  dominant  over  bottom  loss.  This  is  indicated  by.  the  dashed  81  dB 
level  curve  shown  near  the  top  of  Fig.  5.  Rays  with  smaller  N  exhibit  this 
change  in  dominance  at  shorter  ranges.  The  increasing  loss  with  small  range 
values  is  a  consequence  of  any  bottom  reflection  model  which  displays 
monotonicaily  increasing  loss  with  increasing  ray  angle.  For  example,  the 
Rayleigh-3  and  Mackenzie^  models  both  possess  such  a  characteristic. 

The  deviation  from  circles  of  the  level  curves  of  Fig.  5  is  an  influence 
of  bottom  inclination .  For  downslope  propagation  (5  upslope  from  R),  the  ray 
angle  of  incidence  at  the  bottom  becomes  closer  to  grazing  with  each 
successive  bottom  reflection,  the  rapidity  of  approach  to  grazing  depending  on 
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both  a  and  5  bearing  5.  For  8  near  U  and  ±tt  radians,  the  angles  (and 

exhibit  minimum  (and  maximum)  deviation  from  their  horizontal-bottom  values. 

We  observe  that  the  transition  from  bottom  loss  dominance  to  spreading  loss 
dominance  occurs  at  shorter  ranges  for  8  near  0  or  ±tt,  as  indicated  by  the 
closeness  of  contours.  For  example,  there  is  a  10  dB  change  in  transmission 
loss  between  4  and  5  km  range  here,  contrasted  with  a  6  dB  change  when  8  =  ± 
tt/2  rad.  Upslope  propagation  (S  downslope  from  R)  similarly  exhibits  strong 
dependence  on  a  and  8  through  the  bottom-loss  model,  with  the  primary 
distinction  that  the  steepness  of  a  ray  tangent  now  increases  with  each  bottom 
reflection.  This  causes  the  gradient  of  transmission  loss  to  be  smaller  in 
magnitude  when  propagating  upslope  compared  to  downslope.  Thus,  for  a  given 
range,  transmission  loss  for  upslope  propagation  exceeds  that  for  downslope 
propagation.  Qualitatively  similar  level  curves  of  transmission  loss  for 
longer  range  sound  transmissions,  obtained  using  normal  mode  theory,  have 
been  calculated.^ 

The  variation  of  transmission  loss  as  S  assumes  successive  positions  on 
specified  paths  is  shown  in  Fig.  6.  The  paths  P-j  through  P4  appear  in  Fig.  5 
and  are  the  same  as  on  Fig.  2,  with  each  having  a  CPA  of  5  km.  The  dashed 
curve  in  Fig.  6  is  the  loss  for  a  horizontal  bottom  (a=0)  at  depth  dR,  for  a 
linear  path  with  a  CPA  of  5  km.  The  characteristic  dominance  of  bottom  loss 
over  spreading  loss  at  short  ranges  is  easily  seen  in  this  broken  curve.  Path 
P-j  is  an  upslope  track  (downslope  propagation)  mainly  lying  in  a  region 
dominated  by  spreading  loss,  since  bottom  loss  influence  dominates  in  an 
interval  of  about  2  km  on  either  side  of  CPA.  We  note  that  transmission  loss 
on  P1  is  less  than  on  the  same  path  when  a  =  0.  The  downslope  image  of  Pi  is 
P 2?  which  lies  in  a  region  dominated  by  bottom  influence.  The  preeminence  of 
the  efface  of  altitudinal  ray  angle  is  seen  in  the  rapid  increase  of  loss 


as  S  approaches  tt/2  at  CPA.  Overall,,  the  sloping  bottom  causes  an  increase  in 
transmission  loss  and  a  broadness  of  the  region  of  bottom  loss  domination. 

Path  P3  is  directly  upslope,  with  minimum  transmission  loss  occurring  for 
downslope  propagation.  That  part  of  P3  in  which  bottom  loss  dominates,  from 
its  start  to  about  2  km  beyond  CPA,  also  shows  the  strong  influence  of  bottom 
inclination  in  changing  the  spatial  variation  of  transmission  loss  from  the 
dashed  curve.  The  diagonal  trajectory  is  primarily  a  bottom  loss  dominated 
path,  with  spreading  loss  determining  the  quality  of  variation  only  in  the 
last  one  km.  In  summary,  we  see  here  that  a  sloping  bottom  intensifies  the 
variations  in  transmission  loss  for  paths  lying  primarily  in  a  region  of 
bottom  loss  predominance,  as  for  P2  through  P4,  and  reduces  variations  when 
spreading  loss  dominates,  as  for  P-j  . 

We  illustrate  the  influence  of  the  sloping  bottom  by  displaying  relative 
per-rav  transmission  loss,  -20  log-jQ  (A5/A5),  as  level  curves  on  a  spatial 
grid  14  km  square  with  R  at  the  origin  (see  Fig.  7).  The  interval  between 
adjacent  level  curves  is  2  dB  in  this  figure.  When  propagating  downslope  (S 
upslope),  the  bottom  angle  causes  an  increase  in  signal  strength  when  compared 
to  a  horizontal  bottom.  The  converse  is  true  when  signals  are  propagated 
ups lope a  The  strong  influence  of  bottom  inclination  is  seen  from  a  total 
variation  in  transmission  loss  of  more  than  24  d3  over  the  figure.  Regions  of 
small  deviation  from  the  horizontal  bottom  case  are  caused  by  two  different 
mechanisms.  First,  near  the  receiver,  bottom  loss  dominates  because  of  steep 
altitudinal  ray  angles,  and  is  essentially  the  same  in  both  the  unperturbed 
(ct=0)  and  perturbed  cases.  Second,  at  longer  ranges,  there  is  dominance  of 
spreading  loss,  which  has  only  small  deviations,  as  discussed  earlier.  The 


absence  of  an  effect  of  source  bearing  is  seen  near  the  y-axis,  where  little 
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variation  in  0^  occurs  to  produce  changes  in  bottom  loss.  In  the  regions  of 
bottom  loss  dominance,  however,  the  influence  of  6  on  ray  angles  is  expressed 
by  more  closely  spaced  level  curves,  as  seen  for  |  0  |  approaching  tt/2. 

Equations  (9)'  show  that  altitudinal  angles  experience  relatively  large 

change  from  0^,  while  azimuthal  angles  (f)^  vary  only  slightly  from  0,  when 
\&\  is  near  it / 2 .  The  conspicuous  asymmetry  of  the  level  curves  is  caused  by 
bottom  loss  effects  and  the  manner  of  bottom  loss  change  as  discussed  before. 
Thus,  variation  of  0^  with  0  is  the  primary  cause  of  the  large  changes  in 
per-ray  transmission  loss  for  propagation  over  a  sloping  bottom  in  regions  of 
bottom  loss  dominance. 

Before  discussing  total  field,  a  brief  overview  of  the  relationship  among 
its  constituents  is  appropriate.  In  Table  I,  we  display  the  relative 
magnitude  A^/A^  and  the  relative  per-ray  transmission  loss  -20  log^Q  (A^/A^) 
for  the  first  nine  ray  arrivals  at  R  when  S  is  located  upslope  at  (R,0)  = 
(5,-75°)  and  downsiope  at  (5, +75°).  For  these  values  of  0,  5  lies  in  a  region 
where  bottom  loss  dominates  spreading  loss.  As  N  increases,  the  relative 
amplitude  of  an  unperturbed-state  ray  compared  to  that  of  the  strongest  H  =  1 
ray  decreases,  while  -20  log^Q  (A^/AN)  becomes  significant,,  This  shoves  that 
large  changes  in  relative  transmission  loss  are  associated  &dth  ray  components 
of  less  significant  magnitude  in  the  total  acoustic  fieldo  We  see  little 
variation  in  per-ray  transmission  loss  for  N  <*  3;  however*  for  upslope 
propagation  ( S  downsiope ) ,  the  difference  in  bortom  loss  becomes  substantial 
at  N  =  4.  Signals  propagating  downsiope  show  significant  deviation  from  the 
horizontal-bottom  sound  channel  for  N  >_  5.  The  quantitative  differences  with 
N  between  upslope  and  downsiope  propagation  arise  from  the  sensitivity  of  the 
Mackenzie  bottom-loss  model  to  incident  ray  angle  at  the  bottom. 
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We  now  briefly  consider  bottom  slope  effect  on  incoherent  total-field 
amplitude.  Specifically,  we  examine  the  range  R  to  level  curves  of  trans¬ 
mission  loss,  -10  log-jQ  A 2,  where 

A2  =  l  A2  ,  (13a) 

N 

in  which  AN  is  given  by  Eq.  (12).  Similarly,  we  consider  the  constant  range  R 
to  level  curves  of  -10  log-jQ  A^/Where 

A2  =  l  A2  (13b) 

N 

is  the  corresponding  amplitude -squared  for  a  horizontal  bottom  of  depth  dR.  We 
use  the  same  parameter  values  and  S-  R  placement  as  before,  and  consider 
level  curve  values  of  65,  68,  and  71  dB.  For  each  fixed  transmission  loss 
value,  we  calculated  the  percentage  change  in  range,  100  (R-R)/R.  The  results 
appear  in  Fig.  8  as  a  function  of  bearing  angle  8.  We  observe  that  the  range 
to  each  level  curve  is  greater  for  downslope  propagation  (8  <  0)  and  shorter 
for  ups lope  propagation,  so  that  downslope  propagation  exhibits  a  larger 
percentage  variation.  For  the  solid  curve  (71  dB)  ,  R  has  the  constant  value 
6.8  km,  while  R  =  3o8  km  for  63  dB,  and  P  =  1.7  km  for  65  dB.  Percent  changes 
in  range  varies  a  total  of  more  than  40%  as  8  assumes  all  possible  for 

the  71  dB  contour,  reflecting  the  influence  of  bottom  inclination. 

The  primary  cause  of  this  variation  is  the  change  in  ray  angles  at  bottom 
reflection  points  with,  for  downslope  propagation  (i.e.,  8  <  0)  ,  the  rays 
starting  more  steeply  than  for  a  horizontal  bottom  and  becoming  shallower  with 
each  bottom  bounce.  Table  I  shews  that  strong  rays  have  small  N,  whereas  Eqs . 
(9)  imply  that  the  strong  rays  have  least  change  in  geometry  from  the 
horizontal  bottom  rays.  Consequen tly ,  strong  rays  exhibit  transmission  loss 
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close  to  unperturbed  rays.  For  both  the  71  dB  and  68  dB  contours,  the  rays  for 
N  =  1,2,  and  3  dominate  the  total  field;  however,  other  rays  contribute  at 
least  through  N  =  6.  The  lessening  of  the  altitudinal  angle  at  each  successive 
bottom  reflection  causes  range  to  the  fixed  transmission  loss  curve  to 
increase.  Percent  range  variation  for  upslope  propagation  (i.e.,  6  >  0)  is 
explained  similarly.  Rays  with  multiple  bottom  bounces  are  steeper  at 
successive  bottom  reflections.  This  increased  loss  tends  to  shorten  the  range 
to  level  curves  of  transmission  loss.  At  short  range  (i.e.,  the  65  dB 
contour)  ,  both  the  rays  over  a  sloping  bottom  and  those  over  a  horizontal  one 
are  so  steep  that  there  is  little  variation  in  bottom  loss  between  them.  The 
strong  effect  of  S  seen  in  each  curve  is  predicted  by  its  effect  on  each  ray 
from  Eqs .  ( 9 ) . 

III.  COMBINED  SLOPE  AND  FRONTAL  EFFECTS 

Oceanic  fronts  can  occur  in  the  continental  margin  and  in  other  shallow- 
water  regions.30''31  Such  fronts  have  been  modelled  numerically  (see,  for 
example,  Ref.  4)  and  studied  analytically  for  a  horizontal  bottom. 25  ^ 

sloping  bottom  induces  significant  variation  in  acoustical  quantities  as  shown 
in  Sec.  II,  while  Ref.  25  demonstrates  frontally  induced  variations  in  a 
horizontal  bottom  sound  channel.  We  examine  here  the  influence  of  a  plane 
bottom  which  slopes  away  from  shore  on  sound  transmission  through  a  front  in 
shallow  water. 

We  use  a  model  which  defines  the  front  as  a  vertical  plane  separating 
water  masses  of  constant,  but  different,  sound  speed  on  either  of  its  sides. 

Our  simple  model  includes  uniform  but  different  horizontal  currents  on  both 
sides  of  the  front,  as  well.  The  front  is  oriented  parallel  to  the  shore, 
and  we  place  and  R  on  opposite  sides  of  it.  As  in  previous  sections,  the 
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sound  speed  at  R  is  c ,  while  the  souree-side  sound  speed  cs  is  related  to  c 
by 

cs  =  c  +  Ac-  (14a) 

It  is  convenient  to  introduce  the  dimensionless  sound  speed  jump  A,  defined 
by 


A  =  A c/c  .  (14b) 

Typically,  the  along  front  currents  in  shallow  water  vary  in  magnitude  between 
2  and  16  crns""'  ,  but  magnitudes  exceeding  30  cms"^  are  not  uncommon. 31  The 
coordinate  system  established  in  Sec.  I  and  depicted  in  Fig.  1  is  maintained 
here,  and  the  equation  for  the  frontal  plane  is  x  =  Xp.  The  distance  dp  from 
R  to  the  front,  along  the  source-receiver  line,  is 

dp  =  -Xp  esc  0  ,  (14c) 

where  |p|  cannot  be  zero  or  tt  since  we  permit  neither  S  nor  R  to  lie  in  the 
frontal  plane.  The  inclination  angle  0*  of  the  frontal  plane  to  the  S-R  line 
is  related  to  the  source  bearing  by 

0*  =  H(0)ir  -  0  .  ( 14d) 

The  quantities  xp,  dp,  and  0*  are  shown  in  the  inset  of  Fig.  9. 

We  extend  the  notation  of  Sec.  I  to  indicate  the  side  of  the  rent 
(source  or  receiver)  on  which  a  quantity  lies,  by  appending  an  S  or  R  sub¬ 
script.  For  example,  if  the  front  intersects  a  ray  path  between  the  last 
bottom  reflection  and  R,  the  altitudinal  angle  is  G^jg  before  reaching  the 
front,  and  G^nr  a^ter  the  ra-Y  passes  through  the  front.  However,  the 
altitudinal  angle  GRN  of  the  ray  tangent  incident  at  R  is  an  example  of  a 
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quantity  whose  location  is  assumed  known  (on  the  receiver  side  of  the  front), 
so  that  no  additional  subscript  is  written. 

We  determine  ray  geometry  for  fixed  5  and  R  by  tracing  a  ray  forward  from 
S  through  the  front  to  R,  assuming  that  neither  a  surface  nor  bottom  reflec¬ 
tion  point  occurs  at  the  plane  of  the  front.  For  any  given  ray  with  N  bottom 
reflections,  we  take  the  front  to  lie  in  lobe  n  of  that  ray  path.  Ray 
tangents  on  either  side  of  the  front  are  as  described  in  Sec.  I: 


A 

~kNS 


/  cos2ka  cos  0ONS  s^n  §oUS~s^n  ®oNS  sin2ka  \ 
-cos  0ONS  cos  ^oNS 

y  sin  0oNS  cos2ka  +  cos  0oNS  sin  sin2ka  / 
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r  0  ^  k  <c  n  , 
(15a) 
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Two  range  relations  are 


n 

(R-dF)  sin  0  =  £  LkNS  $kNS  (16a) 

0 


and 


N 

dF  sin  0  =  l  LkNR  sin  9  (16b) 

n 

where  is  the  projection  on  the  xy-plane  of  that  portion  of  the  ray  path 

in  lobe  n  on  the  source  side  of  the  front,  and  similarly  for  LhNR  on  the 
receiver  side  of  the  front. 

The  effect  of  the  front  on  a  ray  is  to  bend  it,  the  refraction  relation 
oeinq  — 
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(1+A)An^R  “  ^nNS  “  ^Mp  /  (17a) 

where 

r  =  [(1+A)2  -  i  +  (AnNS*MF)2]1/2  -  AnNS*MF  ,  (17b) 

and  Mp  is  the  unit  normal  to  the  plane  of  the  front  on  the  receiver  side.  As 
in  Sec.  I,  we  can  relate  horizontal  distances,  water  depths,  and  angles  by 
using  Eqs .  (6)  on  either  side  of  the  front.  For  example,  Eq.  (6b)  when 
applied  on  the  source  side  of  the  front  gives  for  1  <_  k  _<  n  as  a  function 

of  h^S'  ®kNS'  ^kNS'  an^  a-  We  insure  continuity  of  the  ray  path  by  matching 
the  ray  traced  backward  from  R  to  the  ray  traced  forward  from  S  using  the 
refraction  relation  Eq.  (17a).  To  trace  from  R,  we  need  one  additional 
relation, 

hNNR  “  (dR  tan  ®NNR  +  °RhR  tan  a  sin  'J’NNr) 

x  ( tan  Qjjnr  -  tan  a  sin  •  (18) 

The  equations  derived  from  Eqs.  (6)  and  Eqs.  (15)— ( 18)  form  a  non-linear 
deterministic  system  for  0oNS,  QNNF'  4>oNS'  and  $NNR* 

Next,  we  write 

®kNS  =  ®N  +  ekNS  >  0  £  k  £  n  ,  (19a) 

®kNR  =  ®N  +  ekNR  •  n  _<_  k  £  N.  ,  (19b) 

4>icNS  =  B  +  XkNS  /  0  £  k  £  n  ,  (19c) 

and 

‘i’kNR  =  3  +  XkMR  '  n  £  k  £  N  ,  (19d) 
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where  the  overbar  on  Q  is  an  altitudinal  angle  of  the  corresponding  ray  when  no 
front  and  no  bottom  slope  are  present.  Thus,  the  e  and  x  terms  represent 
corrections  to  altitudinal  and  azimuthal  angles  which  arise  from  a  sloping 
bottom  and/or  the  presence  of  a  front.  We  assume  that  the  angle  changes  from 
0N  and  8  are  small  relative  to  unity  so  that  we  can  approximate  the  deviations 
by  perturbation  methods.  Further,  the  orders  of  magnitude  of  the  sloping- 
bottom  parameter  a  and  the  dimensionless  frontal  parameters  A  and  signed  Mach 
numbers  =  -  \  vs  |  /c  and  Ur  =  ±  |vR|/c  must  be  specified,  so  that  terms  of 
the  correct  degree  are  kept  in  expansions.  In  particular,  for  moderate  fronts, 
A  =  O(10~2),3°  and  for  a  between  1/2°  and  1°,  a  in  radians  satisfies  A  =  9(a). 
However,  |us|  and  |yR|  are  typically  0(10“4)31  so  that  Mach  numbers  may  be 
ignored  when  only  first  degree  terms  in  A  and  a  are  kept  in  our  expansions. 

Under  these  conditions,  we  approximate  A^s  and  A^nr  from  Eqs.  (15). 
Substituting  the  results  into  Eqs .  (17)  gives  an  approximation  for  F,  and  a 
relation  between  AnNR  and  A nNS  bo  linear  terms  in  a  and  A.  These  results  are 
extended  further  by  additional  approximations  from  Eqs.  (19),  assuming  that 
linear  terms  in  a,  A,  eoNS,  e^NR'  XoNS'  and  Xnnr  are  actually  of  the  same 
order.  We  use  the  ratios  of  the  x  and  y  components,  and  of  the  z  and  y 
components,  from  Eq.  (17a),  to  express  angle  changes  in  terms  of  a.  A,  and 
angle  change  at  S  and  R.  Substituting  the  results  into  Eqs.  (6)  and  (18)  and 
then  into  Eqs.  (16),  and  retaining  linear  terms  in  a  and  A,  we  find  approxi¬ 
mations  of  eoNS,  e^NR'  XoNS'  and  Xnnr  in  terms  of  a  and  A.  Under  the  con¬ 
ditions  that  rays  are  neither  very  steep  nor  close  to  grazing  upon  incidence 
to  the  front  (A  sec^0N  csc^g*  <<  1  )  , nor  c]_ose  to  grazing  on  bottom 


incidence,  we  have 
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eoNS  =  -a{ 2N(N+ohR)R  1  sin  3  cos  0N  sin  0N} 

+  A  dFR-1  tan  0N  ,  (20a) 

eNNR  =  -a{ 2N(N+ohR) R-1  sin  3  cos  0N  sin  0N 

-  2N  sin  3}  -  A( 1 -dFR-1 ) tan  0N  ,  (20b) 


XoNS  =  a  2N(N+cRhR)R  1  cos  3 


and 


+  A  dpR  1  sec^  cot  8*  , 

Xnnr  ®  a[  2N('N+aRhR)R"1  cos  6  -  2N  cos  8  tan  0N] 


-  A  ( 1 -dpR-1 ) sec2  0N  cot  8*  , 


(20c) 


(20d) 


where  dp  =  dp/dp  is  a  dimensionless  distance  from  R  to  the  front.  Recalling 
that  0S^  and  0p^j  are  always  on  the  source  and  receiver  sides  of  the  front, 
respectively,  Eqs .  (1)  complete  the  geometric  approximations.  The  results 
given  in  Eqs.  (20)  reduce  to  those  in  Sec.  I  when  A  =  0*land  reduce  to  those 
in  Ref.  25  [for  Mr  =  Us  =  0]  when  a  =  0. 

We  calculate  travel  time  from 

n  N 

T  =  c"1  E  L  sec  9  +  c_1  7  L  sec  0  ,  (21a) 

N  S  o  kNS  kNS  n  kNR  kNR 
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in  which  the  first  (and  second)  term  is  the  travel  time  from  S  to  the  front 
(and  from  the  front  to  R).  Using  the  approximations  of  Eqs .  (20)  in  Eq.  (21a) 
and  expanding  to  linear  terms  in  a  and  A  gives 

Tjq  =  Tjyj  |  1  +  2Na  sin  B  [N+Oplip^R 
X  [r2  +  (2N+ashs+aRhR)2]  1 

-  (1-dFR-1)A}  ,  (21b) 

where  T^  is  given  by  Eq.  (10c)  and  is  the  travel  time  in  the  absence  of  both 
the  front  and  the  sloping  bottom. 

To  illustrate  the  influence  of  a  sloping  bottom  on  propagation  through  a 

front,  we  calculate  T5  with  c  =  1500  ms””1  ,  Ac  =  20  ms"’1  ,  Xp  =  -2  km,  dR  =  300  m, 

and  source  and  receiver  on  the  surface.  The  value  of  Ac  is  typical  of  the  Slope 

Front  in  the  North  Atlantic.  32  The  results  for  travel  time  AT5  =  T5  -  T5  are 

shown  in  Fig.  9,  for  three  values  of  bottom  slope  a,  as  S  occupies  successive 

positions  along  two  paths,  each  with  a  CPA  5  km  downslope  from  R.  The  solid 

path  P<|  is  parallel  to  the  shore,  whereas  the  dashed  path  p?  tracks  diagonally 

upslope.  The  curves  for  a  =  0  represent  the  influence  of  the  front  in  a 

horizontal  bottom  sound  channel.  Note  that  on  path  P-j  .  the  front  induces 

maximum  IAT5I  at  the  beginning  (and  end)  of  P-j  where  3  is  closest  to  tt  (or  zero) 
and  where  the  relative  front  location  dpR  1  is  smallest.  For  a  =  0.5°  the 

bottom  inclination  cancels  the  frontal  effect  at  the  ends  of  the  track  and 

increases  AT5  by  more  than  60  ms  at  CPA  over  that  for  a  =  0° .  When  a  =  1°,  the 

bottom  angle  accounts  for  more  than  140  ms  of  the  value  of  AT5  at  CPA. 
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For  5  tracking  diagonally  upslope  on  path  ,  symmetry  about  CPA 

disappears.  At  the  end  of  the  track  when  a  =  0° ,  AT5  is  20  ms  more  than  its 

value  at  the  start  of  the  path.  For  larger  a  we  again  see  an  increasingly 

strong  influence  oF  the  sloping  bottom-  Notice  that  the  sloping  bottom  exerts  a 

weaker  influence  near  the  end  of  P2  because  0  is  becoming  small,  causing  the 

inclination  induced  travel-time  change  to  be  small,  as  was  suggested  in  Fig.  4. 

The  strongest  frontal  effect  occurs  when  relative  front  location  dpR-1  is 

smallest  at  the  start  of  P2/  but  the  inclined  bottom  effect  is  large  there  also, 

because  0  is  near  tt/2 .  Thus,  we  see  that  the  peak  of  AT5  when  a  *  0  occurs  before 

CPA  because  of  the  asymmetric  variation  in  frontally  induced  travel  time  change. 

A  primary  qualitative  difference  in  the  behavior  of  AT5  on  the  two  paths 

comes  from  the  manner  in  which  the  front  affects  travel  time.  When  a  =  0,  Eq. 

(21b)  predicts  the  behavior  we  illustrate.  Frontally  induced  travel  time  change 
varies  as  d?R“  1  and  shows  symmetric  behavior  about  CPA  on  P^  ,  but  monotonicity  on 

P2-  The  change  in  AT5  induced  by  the  sloping  bottom  alone  on  path  P-j  exhibits  the 

same  symmetric  increase  and  decrease  about  the  CPA  as  did  the  travel-time  change 

caused  by  the  front.  These  effects  are  additive  when  a  £  0  and  produce  curves 

with  the  same  general  properties.  In  contrast,  the  variation  in  AT5  caused  by  the 

sloping  bottom  alone  is  asymmetric  about  the  CPA  and  has  an  increasing  and  then 

decreasing  behavior.  When  a  *  0,  the  asymmetry  is  more  pronounced,  with  the 

dominant  influence  of  the  sloping  bottom  determining  the  illustrated  behavior. 

The  influence  of  the  front  on  AT5  is  more  pronounced  for  a  =  0.5°  than  for 

a  =  1°  because  the  relative  magnitude  of  frontally  induced  AT5  (a=0)  compared  to 

that  caused  by  the  bottom  angle  is  less  when  a  =  1°. 

Downslope  propagation  is  illustrated  by  placing  S  upslope  from  R,  and 

adjusting  the  frontal  parameter  Ac  to  correspond  to  warmer  water  offshore  from  the 

front.  We  display  AT5  for  the  same  three  values  of  a  on  two  paths  upslope  from  R 
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in  Fig.  10,  with  c  =  1520  ms"1,  Ac  =  -20  ms-1,  xF  =  2  km,  and  other  parameters  as 
in  Fig.  9.  The  paths  are  similar  to  those  previously  described,  except  that  is 
parallel  to  the  shore  with  CPA  5  km  upslope  from  R.  Path  P2  again  tracks  diagon¬ 
ally  upslope  with  a  5  km  CPA  to  R.  Since  we  are  now  transmitting  from  cold  to 
warm  water  (Ac  <  0),  the  front  alone  (a  =  0)  causes  AT5  to  be  positive.  However, 
as  seen  in  Fig.  4,  the  sloping  bottom  generates  negative  AT5  for  downslope 
propagation.  The  effect  of  the  sloping  bottom  on  the  frontal  influence  for 
downslope  propagation  is  analogous  to  that  for  upslope  propagation.  Comparison  of 
Figs.  9  and  10  reveals  that  the  concavity  reversal  from  upslope  to  downslope 
propagation  is  a  result  of  Ac  changing  sign,  with  a  corresponding  change  in  the 
sloping  bottom  travel  time  effect.  The  difference  in  locations  of  peak  (minimum) 
AT5  when  a  r  0  on  path  P2  is  a  result  of  the  different  locations  of  peak  (minimum) 
AT5  caused  by  the  sloping  bottom,  as  suggested  in  Fig.  4.  Traversing  path  P2  in 
Fig.  1 0  in  an  opposite  direction  causes  the  relative  extrema  of  AT5  to  occur 
slightly  before  CPA  as  in  Fig.  9.  The  strong  travel  time  influence,  on  signals 
transmitted  through  a  shallow  water  front,  of  the  bottom  inclination 
(approximately  150  ms  shown)  is  qualitatively  similar  whether  transmitting  upslope 
or  downslope. 

As  discussed  in  Sec.  II  and  in  Ref.  24,  the  sloping  bottom  exerts  a  much 
stronger  effect  on  transmission  loss  than  that  produced  by  a  front.  Conse¬ 
quently,  we  will  not  consider  those  combined  effects  here.  Results  of  this 
section  have  shown  that  bottom  inclination  has  a  profound  effect  on  travel 
time  changes  induced  by  a  moderate  strength  front  in  shallow  water. 

IV.  SUMMARY 

Effects  of  a  sloping  ocean  bottom  on  three-dimensional  per-ray  arrivals 
at  a  receiving  point,  and  on  incoherent  total-field  transmission  loss,  are 
examined.  The  ocean  medium  is  usually  modelled  as  isospeed,  although  an  ocean 
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with  two  distinct  sound  speeds  is  used  to  examine  propagation  through  a  front 
on  the  continental  shelf  and  parallel  to  the  shore.  The  sound  source  and 
receiver  are  at  arbitrary  depths,  and  their  relationship  to  the  sloping  bottom 
is  arbitrary.  Short  range  propagation  for  source-receiver  paths  having  both 
upslope  and  downslope  components  are  investigated.  New  approximations  for  ray 
geometry  over  the  sloping  bottom  are  found  by  perturbing  from  the  basic  state 
in  which  the  bottom  is  horizontal  with  water  depth  of  the  receiver. 

Principal  features  of  variations  induced  by  the  sloping  bottom  are 
investigated.  An  approximation  to  change  in  per-ray  travel  time  is  derived 
which  demonstrates  dependence  on  source-receiver  location  over  the  slope  and 
on  bottom  angle.  In  numerical  examples,  it  is  shown  that  the  travel  time 
change  for  a  ray  with  five  bottom  reflections  is  approximately  ±  200  ms  at  a 
range  of  about  5  km,  where  the  sign  depends  upon  source-receiver  orientation. 

One  of  the  new  results  from  our  formulas  is  that  although  the  sloping 
bottom  causes  an  insignificant  change  in  spreading  loss,  it  induces  significant 
changes  in  per-ray  transmission  loss  through  bottom  loss  effects.  This  fact  is 
illustrated  for  a  Mackenzie  bottom,  but  it  would  hold  for  many  other  bottom 
models.  Upslope  propagation  is  shown  to  experience  greater  transmission  loss 
than  downslope  propagation.  In  comparison  to  a  flat  bottomed  sound  channel,  a 
bottom  inclination  angle  of  2°  induces  a  decrease  (or  increase)  of  approxi¬ 
mately  10  dB  (or  14  dB)  when  propagating  downslope  (or  upslope)  o^ier  a  range  of 
about  5  km. 

The  range  from  receiver  to  level  curves  of  total-field  transmission  loss 
is  examined  as  a  percentage  change  from  the  range  to  corresponding  level 
curves  in  a  sound  channel  with  horizontal  bottom.  It  is  found  in  examples 
that  the  range  varies  by  over  ±  20^,  the  sign  depending  upon  down  or  upslope 

i 

propagation . 
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Finally,  we  add  to  the  sloping  channel  a  model  of  a  moderate  strength 
shallow  water  front,  parallel  to  the  shoreline,  and  characterized  by  different 
sound  speeds  on  either  side  of  a  vertical  plane.  Geometric  and  travel  time 
approximations  are  developed  when  bottom  slope  and  relative  change  in  sound 
speed  across  the  front  are  of  the  same  magnitude.  We  find  that  bottom 
inclination  induces  significant  variations  in  the  travel  time  change  predicted 
by  the  frontal  model  with  a  flat  bottom.  Indeed,  numerical  examples  are 
presented  in  which  bottom  slope  induces  a  400%  increase  in  travel  time  change 
when  a  signal  is  propagated  upslope  over  a  range  of  5  km. 
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TABLE  I.  Relative  amplitude  and  transmission  loss  (dB)  for  first  nine 
rays  incident  on  R.  R  =  5  km,  S  =  75°  (  S downslope ) , 

6  -  -75°  (S  upslope);  other  parameters  as  in  Fig.  2. 


N 

An/A1 

Rel.  Magnitude 

-20  log -j  q  ( A^j/A^ ) 

Upslope  Downslope 

1 

1  .00 

-  0.1  1 

+  0.09 

2 

0.91 

-  0.43 

+  0.26 

3 

0.80 

-  0.55 

+  0.94 

4 

0.66 

-  1.15 

+  8.08 

5 

0.31 

-  4.98 

+14.00 

6 

0.042 

-14.96 

+  9.41 

7 

0.009 

-18.77 

+  7.41 

8 

0.002 

-20.39 

+  6.47 

9 

0.0007 

-20.92 

+  5.24 
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